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The fluctuation theorem establishes general relations between transport coefficients and fluctua¬ 
tions in nonequilibrium systems. Recently there was much interest in quantum fluctuation relations 
for electric currents. Since charge carriers also carry spin, it is important to extend the fluctuation 
theorem to spin currents. We use the principle of microscopic reversibility to derive such theorem. 
As a consequence, we obtain a family of relations between transport coefficients and fluctuations of 
spin currents. We apply the relations to the spin Seebeck effect and rectification of spin currents. 
Our relations do not depend on a microscopic model and hence can be used to test the validity of 
theoretical approximations in spin-transport problems. 
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I. INTRODUCTION 

Both charge and spin of electron play major roles in 
solid state physics and its applications but the electric 
and spin currents behave in different ways. Coulomb 
interaction brings a high energy cost to any violation 
of charge neutrality. This fact and charge conservation 
make it easy to create and manipulate charge currents. 
At the same time, spin does not conserve. Its coupling to 
external probes is weaker than that of the electric charge. 
Thus, manipulating spin currents is a harder task than 
working with electric currents. Recently, much effort has 
focused on overcoming this challenge in the hope to get 
more control over the spin degree of freedoroi. That ef¬ 
fort has been motivated by both potential applications 
and the interest in basic spin physics. An important 
open problem in the field is the understanding of non¬ 
linear spin transport. 

The problem is challenging because nonlinear trans¬ 
port occurs away from thermal equilibrium. Until re¬ 
cently, any progress in the theory of nonequilibrium 
transport was inevitably based on microscopic models^. 
The situation changed after the discovery of fluctuation 
theorems^^— which brought a powerful general principle 
to nonequilibrium statistical mechanics. The goal of this 
paper consists in the derivation of a family of fluctu¬ 
ation relations for spin currents and noises. We then 
apply them to several spin-transport problems, includ¬ 
ing rectification of spin currents and the spin Seebeck 
effecliSi^— , i.e., the generation of a spin current by a tem¬ 
perature gradient. An important application is the the¬ 
ory of spintronic nanomachines^^. Our fluctuation rela¬ 
tions are completely general and do not depend on mi¬ 
croscopic details or the choice of a model. This makes 
their application particularly interesting. They also give 
a useful tool for testing theoretical approximations. 

The history of fluctuation relations is fascinating. 
They were first introduced by Bochkov and Kuzovlev in 
the late 1970s^“— . The name “fluctuation theorem” was 
coined by Hanggi in 1978^^. Yet, its importance has been 


fully appreciated only after a series of seminal results in 
the 1990s, in particular, the works by Evans, Cohen and 
Morrissi^, Jarzinskii^, and Crooksi^. Fluctuation the¬ 
orems generalize the principle of detailed balance. De¬ 
tailed balance establishes a fundamental relation between 
the rates of different processes. Fluctuation theorems go 
a step further and relate the probabilities of a process 
and its time-reversed version on a macroscopic time scale. 
At first sight, fluctuation relations apply to time-reversal 
invariant systems only. Indeed, a process and its time- 
reversed version can occur in the same system only if the 
time-reversal symmetry is present, i.e., in the absence 
of a spontaneous magnetization and external magnetic 
fields. Otherwise, the time-reversed process would re¬ 
quire changing the sign of the magnetic field. However, 
interesting relations do not have to connect the proper¬ 
ties of a single system. This fact was already recognized 
in the pioneering papers^i^. A number of fluctuation re¬ 
lations for electric currents in a magnetic field were de¬ 
rived in Ref. [O- Those relations contain currents and 
noises in two systems that differ by their directions of 
the magnetic fields. In some cases it even proved possi¬ 
ble to derive fluctuation relations that connect currents 
and noises in a single system in the absence of the time- 
reversal symmetryi^d^. A review of those developments 
can be found in Ref. [^O- 

Focusing on the spin degree of freedom means treating 
currents as multi-component with the components, cor¬ 
responding to the up and down spin projections. A fluc¬ 
tuation theorem for multi-component currents was ad¬ 
dressed in Ref. [2ll That theorem is not directly relevant 
in our problem since Ref. considers the currents of 
time-reversal-invariant quantities while spin changes its 
sign under time reversal. More specifically, the spin op¬ 
erator S satisfies 


0S0-1 = -S, (1) 

where 0 is the time reversal operator. In contrast, charge 
is invariant under time reversal: 0(50“^ = Q. Several 
fluctuation relations for spin currents were derived in Ref. 
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l22l . In this paper we establish a bigger family of relations 
which could not be obtained with the approach of Ref. 
[^ . Indeed, Ref. did not make use of microreversibility 
that allows us to go beyond the relations^^. 

According to the principle of microscopic reversibility, 
the amplitudes of a process in a given system and the 
time-reversed process in the system with the reversed 
directions of all magnetic fields must be complex conju¬ 
gate. On the basis of calculations within the Landauer- 
Biittiker formalis m^^i^'^ , it was proposed that microre¬ 
versibility fails in nonzero magnetic fields. We believe 
that the conflict of the Landauer-Biittiker approxima¬ 
tion and the principle of microscopic reversibility shows 
only the limitations of the Landauer-Biittiker approach. 
Note that the conclusions, based on the microreversibil¬ 
ity principle are supported by the experimen t^^'^^ and 
agree with theoretical calculations beyond the Landauer- 
Biittiker approximation2Zi2^. A detailed discussion of mi¬ 
croreversibility in the presence of a magnetic field can be 
found in the review^S. To make this paper self-contained 
we repeat some of that discussion in AppendixjA] Clearly, 
it is important to address the consequences of microre¬ 
versibility for spin currents and noises. 

Besides Refs. several other articles address 

related physics. In particular. Refs. explore fluc¬ 

tuation relations for time-dependent protocols that can 
break time-reversal symmetry. In contrast to our work 
and similar to Ref. they do not consider observables 
that change their sign under time-reversal. 

The connection of our results with the pioneering work 
is particularly interesting since these works consider 
quantities, even and odd under time-reversal, on equal 
footing. Refs, yja contain great wealth of results. In 
particular, Ref. Isj derives general fluctuation relations 
in thermal equilibrium in the absence of magnetic fields. 
Our results are related to those of Ref. H in the equi¬ 
librium limit at zero magnetic field. At the same time, 
our main focus is on nonequilibrium transport, driven by 
the differences of the reservoir temperatures and chemical 
potentials, in an external magnetic field. 

Refs. 00 derive fluctuation relations for nonequilib¬ 
rium steady states that form after a sudden change of 
parameters in a system which is initially in equilibrium 
and does not break time-reversal symmetry (see Section 
7 of Ref. i). This problem is related to but different from 
ours. Indeed, the meaning of the fluctuation relations^ is 
different from our work. Ref. 0 allows one to connect the 
derivatives of various correlation functions with respect 
to the parameters of the Hamiltonian. This is relevant in 
the problem of quantum quench (see, e.g., Ref. M)- At 
the same time, our focus is on the dependence of corre¬ 
lation and response functions on the reservoir tempera¬ 
tures and chemical potentials. This is a natural question 
in the context of mesoscopic transport^^. Note that the 
approach^ has been used to derive fluctuation relations 
for electric currents, driven by a voltage bias between the 
reservoirs of the same temperature. This proved possible 
because electro-neutrality dictates that a change of the 


potential energy of a reservoir leads to an equal change 
of its electro-chemical potential. It is unclear if a similar 
approach can be used in our problem which involves gra¬ 
dients of the temperature and of the chemical potential, 
conjugate to the spin. We use a different method that 
builds on Refs. 1^1^ . 

Our paper has the following structure. We set the 
problem and derive a general fluctuation theorem for spin 
transport in Section IE In the third Section we deduce 
from the theorem a family of fluctuation relations for spin 
currents and noises. Section nYi applies our relations to 
the problems of spin-current rectification and to multi¬ 
terminal spintronic devices. We compare our results with 
the previous work in subsection IIVDI Section |V] focuses 
on thermospin transport and Section IVTl summarizes our 
results. In Appendix lAl we discuss the microreversibility 
controversy. Appendix [E reviews ways to measure spin 
currents and chemical potentials. Technical details of 
our calculations are relegated to Appendices O and |E 
Appendix [E] addresses spin transport in chiral system, 
where stronger and simpler fluctuation relations can be 
derived. 


II. FLUCTUATION THEOREM 

We will study fluctuation theorems that are concerned 
with the transport through a finite system (conductor) 
attached to two or more large reservoirs (Fig. [T|). Each 
reservoir is assumed to be in thermal equilibrium but the 
reservoirs are not necessarily in equilibrium with each 
other. A difference between the reservoirs’ chemical po¬ 
tentials drives a current through the conductor. Our 
goal is to establish various relations between the correla¬ 
tion and response functions of the currents. The crucial 
assumption is that the system does not have long-term 
memory beyond a characteristic transport time through 
the conductor. Thus, we can assume any convenient 
choice of the system’s past. The standard choice is as 
follows. Initially, the reservoirs and the conductor are 
disconnected. Next, the particle and energy exchange 
between the system and reservoirs is turned on for some 



FIG. 1: Sketch of the A-terminal setup (only four terminals 
are shown) for spin, charge and heat transport. 
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long time r. At the end, the reservoirs are again dis¬ 
connected from the conductor. We are interested in the 
transport during the time interval r. Yet, it is most con¬ 
venient to describe such transport in terms of the quanti¬ 
ties, measured before the interaction between the conduc¬ 
tor and the reservoirs was turned on at t = 0 and after 
the interaction was turned off at t = r. For example, 
the average current from reservoir number k is dehned as 
— [Qk{T) — (5fc(0)]/T, where Qk{t) is the total charge of 
reservoir k at time t. The energy and spin currents are 
defined in a similar way. Thus, the following measure¬ 
ment protocol is assumed: we measure the charge, the 
z-projection of the spin and the energy of each reservoir 
before and after the evolution period r. 

Such simultaneous measurements are only possible, if 
the charge and the z-component of the spin commute 
with the energy operator and hence conserve. This poses 
an obvious problem: spin does not conserve because of 
the spin-orbit interaction. In addition, external mag¬ 
netic fields may lead to the Larmor precession which 
also breaks spin conservation. Moreover, if the spin- 
relaxation time is short it may be meaningless to even 
speak about spin currents, flowing from the reservoirs. 
Indeed, the spin current is defined as the total spin flow¬ 
ing per unit time through a cross-section near the bound¬ 
ary of the reservoir and conductor. In the absence of spin 
conservation such current depends on the choice of the 
cross-section. 

The above discussion necessitates the following as¬ 
sumptions. First, to avoid the problem with the Lar¬ 
mor precession, we assume that the magnetic fields are 
oriented along the z—axis in the reservoirs. Such mag¬ 
netic fields would break the conservation of the x- and 
y-components of the reservoir spin but we are only in¬ 
terested in the z-component. Second, we assume a weak 
spin-orbit interaction in the reservoirs. The interaction 
must be weak enough that it can be neglected on the time 
scale of the transport through the conductor and on the 
much longer time scale during which a quasiequilibrium 
establishes in a large portion of each reservoir near the 
interface with the conductor. Such quasiequilibrium can 
be described by two different chemical potentials in each 
reservoir for the up and down projections of the elec¬ 
tron spin on the z-axis. In the simplest situation the two 
chemical potentials are equal but we will consider the 
most general case of arbitrary chemical potentials. We 
will assume that the quasiequilibrium survives during the 
evolution period r. Note that very long spin-relaxation 
times were reported for donor spins in pure silicon^. 

The reservoir spin is gradually lost due to the spin- 
orbit interaction and the spin current into the conduc¬ 
tor. Hence, the two chemical potentials slowly depend 
on time. To keep them constant, one would need to com¬ 
pensate the spin loss by a spin current injection into the 
reservoir. 

We do not make any assumptions about the spin-orbit 
interaction and the direction of the magnetic field in the 
conductor. If the z-projection of the spin does not con¬ 


serve in the conductor then the spin currents, leaving the 
reservoirs, do not add up to zero. 

We are now ready to derive the fluctuation theorem. 
Our method is connected with the approach of Refs, 
mil. We assume that initially all reservoirs are discon¬ 
nected. It will be convenient to include the conductor as 
a part of one of the reservoirs. This requires us to make 
the Hamiltonian of the conductor time-dependent so that 
the spin conserves in the conductor at t < 0 and t > t. 
We will assume that the time-dependence is slow (i.e., 
T is large) and hence the energy conserves. We will use 
the units such that ft = fcs = 1. We also assume in this 
section that the temperature is the same in all reservoirs. 
Spin currents, driven by temperature gradients, will be 
addressed in Section V. 

Specihcally, we assume the following protocol. The 
system is described by the time-dependent Hamiltonian 
H = where N is the number of 

the reservoirs. Hi are time-independent Hamiltonians of 
the reservoirs, ifto(t) is the Hamiltonian of the conductor, 
and H[)i (t) describes the interaction of the conductor with 
reservoir i. We asuume Hoi{t) = 0 at t < 0 and t > 
T. Also, HQi(t) are time-independent at times t such 
that tq < t < T — tq, where tq <g; r is the interaction¬ 
switching time. To should be large enough so that we 
can employ energy conservation. We also assume that 
Hq does not depend on time at t < 0 and t > t, and 
Ho{t < 0) = Hoit > t). Its precise form is unimportant 
at such t but it is assumed that at t < 0 and t > t the 
z-component of the spin conserves in the conductor, i.e., 
the spin-orbit interaction is absent. At tq < t < r — tq, 
Hq is time-independent and equals the experimentally 
relevant Hamiltonian of the conductor. In what follows 
we absorb all parts of the system into various reservoirs 
and thus redefine the reservoir Hamitonians as Hi(t) = 
Hi + Hoi{t) + Ho{t) and Hk{t) = Hk + Hok{t), k > 1. 
Each reservoir is initially {t = 0) in equilibrium. Different 
reservoirs are not in equilibrium with each other, i.e., 
their chemical potentials are not the same. Since the 
conductor is absorbed into the first reservoir, its initial 
state is equilibrium with the same chemical potentials 
as in the first reservoir. Certainly, the finite size of the 
conductor means that the choices of its initial state and 
initial Hamiltonian do not affect our results. 

The initial density matrix of the system is the product 
of the density matrices of the reservoirs: 


Pf = Y[p^ = 

i 


n 


1 

^exp 


V^Q^ + Vs^S.^ - E, 

T 


( 2 ) 


where 


Zi = Tr exp 


V^Q^ + Vs^S,^ - E, 

T 


( 3 ) 


Qi is the charge of reservoir i, Ei is its energy, S^i is the 
z-component of the spin in reservoir i, Vi is the electro¬ 
chemical potential and Vsi is the spin chemical potential. 
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conjugate to Szi (experimental procedures for the detec¬ 
tion of spin currents and spin chemical potentials are dis¬ 
cussed in Appendix IbI) . The chemical potentials for the 
up and down spins are {—e)Vi ± Vsil2 and equal —eVi 
at Vsi = 0, where —e is the electron charge. Eq. © 
assumes that all reservoirs have the same temperature 
T. 

The evolution operator is 

UF{t) = T exp f—i f (4) 


where H{t) = ^ Hi (t) is the time-dependent Hamilto¬ 
nian, and T is the time ordering operator. 


We can now write the probability of the forward pro¬ 
cess. This is the probability that during the evolution 
period r the charges of the reservoirs change by AQi 
and their spins change by ASzi- 


PF{AQ,ASz)=Y.{Q{0),Sz{0),n{0)\pF\Q{0),Sz{0),n{0)) x |(Q(r),S,(r),n(T)|[/;^(r)|Q(0),S,(0),n(0))|' 

X <5(AQ - (Q(r) - Q(0)))<5(AS, - (S,(t) - S,(0))), (5) 


where A stays for a vector whose component Ai describes 
reservoir i. The summation in Eq. extends over all 
initial and final states which we label with the charge 
vector Q, the spin vector and an additional quantum 
number n that distinguishes quantum states with identi¬ 
cal charges and spins in the reservoirs. 

We now introduce the time-reversed backward process. 
The initial density matrix 


Pb = 0pF0 ^ 


n 


1 

:^exp 



Vs.Sz^ - E, 

T 


( 6 ) 

where 0 is the antiunitary time-reversal operator, satisfy¬ 
ing 01 = — 10 , and we use the relation QSziQ~^ = —Szi- 
The dynamics is controlled by the Hamiltonian 


Hb =QH{t (7) 


The evolution operator 


exp i J HB{t')dt‘ 

=0 T“^exp(-|-i f H{t')dt') 

Jr — t 


0 - 


( 8 ) 


where T ^ is the reversed time ordering operator. Hence, 

[/b(t) = 0t/;(r)0-\ (9) 

This equation expresses the microreversibility principle. 

We compute the probability that during the evolution 
period r the changes of the reservoir charges are given 
by the vector —AQ and the reservoir spins change by 
+AS,: 


Pb(-AQ, AS^) = ^[(Q(T),S^(T),n(T)|0 ^]pB[0|Q('r),S2(T),n(T))] 

X |[(Q(O),S,(O),n(O)|0-i]C/B(T)[0|Q(r),S,(T),n(T))]|'5(-AQ - (Q(0) - Q(r)))<5(AS, - (S,(0) - S,(r))). (10) 

The above formula differs from Eq. ([5]) in that r stays in place of 0 and vice versa inside the bra and ket vectors. 
This change of notation emphasizes the time-reversed dynamics. For the same reason, all matrix elements are written 
in the basis of the time-reversed states 0|Q,S;;,n). Next, we combine Eqs. (O, (flUl) . (I2|), (O and dH]) and use the 
conservation of the energy ^ Ei on the evolution period r. We find 


Pf(AQ, as,) = Pb(-AQ, as,) exp 


'^[ViAQi + VsiASzi) 
T 


( 11 ) 


Finally, we note that Pb can be interpreted as the probability of a forward process in the system with the same 
chemical potentials V and the opposite chemical potentials — Vg in the opposite magnetic field. Here, the fact that 
V 5 changes its sign follows from the comparison of Eqs. m and ®. Strictly speaking, in addition to the sign 
changes of the magnetic field and V 5 , one also needs to demand the opposite sign of the reservoir magnetizations in a 
system with spontaneously broken time-reversal symmetry. In what follows we instead assume that the time-reversal 
symmetry is always broken explicitly by some (possibly infinitesimal) magnetic field. Interpreting Eq. dill) as a 
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relation between the probabilities of the two forward processes and getting rid of the indices B and F, one obtains 
the fluctuation theorem 


P(AQ, AS,, V, Vs, B) = P(-AQ, AS,, V, -Vs, -B) exp 


J:{V,AQ, + Vs^AS,,) 


T 


( 12 ) 


where B is the magnetic field. The conservation of S', in 
the reservoirs implies that the magnetic field is uniform 
inside each reservoir. Otherwise, the Maxwell equations 
would require the field to have a nonzero x ov y com¬ 
ponent. The field does not have to be uniform in the 
conductor. 

It what follows it will be convenient to slightly modify 
the notations. We will redefine Vi ^ Vi — V, where V 
is the common electro-chemical potential of all reservoirs 
in equilibrium. This will not change the form of Eq. 
(fT^ but will allow us to assume that V are small in the 
limit of the low voltage bias. On the other hand, Vsi 
are already small in the limit of the low voltage bias since 
spin is not conserved far from the conductor and hence its 
chemical potential is zero in a true thermal equilibrium. 


III. FLUCTUATION RELATIONS FOR 
TRANSPORT COEFFICIENTS 


The fluctuation theorem (fT^ is a relation for the prob¬ 
ability distribution functions which are hard to extract 
from experiment. Transport data, such as electric and 
spin currents, conductances and noises are among the 
most important observables. In this section, we estab¬ 
lish several relations {fluctuation relations) among these 
transport quantities from the fluctuation theorem (II2D . 


A. Definitions of transport quantities 

Let us begin with the definitions of the transport pa¬ 
rameters. Recall that our system remains in a steady 
state for most of the protocol duration r (Sec. mi). Un¬ 
der this assumption, the average electric current fl and 
the average spin current Isi, injected into the conductor 
from reservoir i, can be written as 

h = — lim — (A(3i)B,v.Vs, (13) 

T—¥00 T 

Isi = — lim — (AS',i)B,v.Vs, (14) 

T—^■OO T 


where the brackets (... )b.v,Vs mean averaging over the 
distribution function P(AQ,AS,,V,V 5 ,B). When no 
ambiguity arises, we will simply use (...) instead of 
(• • • )b,v,Vs for brevity. In the limit of low voltages, one 


can expand and Isi as 


+ ^Vs, 


dVs, 


disr. 




b ^3 ^ QY^ 3 


2^dVdVk ^ 

jk 

( 15 ) 


2 ^ dVsjdVsk 

JK 


(16) 


where the derivatives are taken at V = V5 = 0. The co¬ 
efficients Gij and are the usual multi-terminal con¬ 
ductances and spin conductances. The coefficient 
is the linear response of the electric current li with re¬ 
spect to a small change in Vsj, and similarly for 
The higher-order derivatives are nonlinear response coef¬ 
ficients. 

We are also interested in various cumulants (correla¬ 
tion functions) of the currents. In particular, the noise 
of the electric current is defined as 

= 2 J dtdt'ihit) - h){Iflt') - I,)/t 
= 2 J dtdt'[Ii{t)Ij{t') - Iilj]/T, (17) 


where the bar means the average with respect to quantum 
and thermal fluctuations, and the factor of 2 follows a 
convention. The noise of the spin current and the 
cross-noise between the spin and electric currents are 

= 2 J dtdt'lJ^i)!^ - Js^isj]/T, (18) 

Sflj = 2 J dtdt'[F{t)Isj{t') - UsA/t. (19) 

The third cumulant of the electric current 

a, AJ 

= — J dtdt'dt" ^liljlk — [flljlk + Ijlkh + Ikfllj] 

+ 2lJjIk). ( 20 ) 

We would like to emphasize that the third cumulant is 
accessible experimentally^^. 

It is convenient to rewrite the above definitions in 
terms of AQi and AS',^. The cross-noise between fl and 
Ij can be written as 

S,, = lim - ((Aft - (AQ,))(AQ, - (AQ,))). (21) 

T—¥00 T 
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Of particular interest are lower-order coefficients in the 
Taylor expansion of in the powers of V and V 5 . 
Keeping up to the linear terms in V and Vs, we can 
approximate Sij by 

Sij = lim -{AQiAQj) ( 22 ) 

T—¥00 T 

To obtain (1^ one uses the fact that the average currents 
vanish at zero bias, i.e., (AQi)/r is linear in V, Vs in the 
limit of small V, Vs and large r. 

Similarly, in the linear order in V,Vs, the noises of 
the spin currents Isi and Isj and the cross-noises of li 
and Isj can be expressed as 

Sf^ = lim -{AS,,AS,,), (23) 

S^] = lim -{AQ,AS„). (24) 

In the zeroth order in the bias, the third cumulants of 
the spin and electric currents can be written as 

Cijk = - lim -{AQ,AQjAQk), (25) 

T —>00 7" 

= - lim -{AS„AS„AS,k). (26) 


In the bulk of this paper we will focus on the relations 
between the currents Ea. (ll3ll4l) . noises (I22II24I) in the lin¬ 
ear order in the bias, and the third cumulants (1251 
at zero bias. One can also define higher-order cumulants 
and keep their dependence on V, Vs to an arbitrary or¬ 
der. However, they are experimentally less relevant. We 
discuss higher-order cumulants in Appendix [C] for com¬ 
pleteness. 


B. Derivation of fluctuation relations 


We are now ready to use the fluctuation theorem (fT^ 
to establish several fluctuation relations for the transport 
coefficients. These relations include the famous Nyquist 
formula and the Onsager reciprocal relations for electric 
currents, and their generalizations to spin currents. 


I 

To begin, we introduce the following piece of notation: 


(B,V,Vs)= lim -(TTAQ^'^A^^Ob.v.v. = lim 

r—^00 T T—^oc 


- E 

T—^■OO T 

AQ,AS* 


n^(AQ, AS,, V, Vs, B). (27) 


This notation allows a compact and simple formulation of the expressions (fT^ [T4l [221IM1) . Higher-order cumulants 
are expressed in terms of combinations of the above quantities (127(1 with different n, m. 

We next consider the Taylor expansion of C"'™(B, V, Vs) near V = Vs = 0: 

C"’-(B,V,Vs) = ^TD;;r(B)n(^) ' ■ ( 28 ) 


The coefficients L";™(B) are the transport coefficients that we are interested in. For example, the conductance G^- 
equals —rL"’™(B) with n, = i/j = I and all other entries of n,m, iz,/x vanishing. The spin conductance G^ equals 
—rL"^™(B) with rrii = = 1 and all other entries vanishing. The physical meaning of the other coefficients L"’*” 

is seen from their relation with various noise cumulants. For example, Fq. (l26l) . where the three 

nonzero entries in niijk are rrii = rrij = mk = I. 

Next, we substitute the fluctuation theorem (I12p in the right hand side of Fq. (iTTi) and Taylor expand the right 
hand side in powers of ^ and A comparison of this Taylor series with the expansion (1281) yields a family of 
general relations among the transport coefficients at the opposite orientations of the magnetic field 


i":;r(B) = (-i)^> 


{rii+fii) 


EEn 


U,\Wi\ 


j n+u,m-[ 
u,/x— 


'(-B) 


(29) 


where u — u denotes the vector with the components Vi — Ui, and the summation X(u=o S«(=oSu 2 =o '' ’ • 
the following sections, we will use Fq. (1291) to extract fluctuation relations for experimental observables in various 
settings. 


IV. EXAMPLES the experimentally most relevant transport quantities, in- 

We now apply the fluctuation relations (EH) to sev¬ 
eral physical systems. Our focus is on relations between 
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eluding currents, noises and the third cumulants. In this 
section we assume that all reservoir temperatures are the 
same. The spin Seebeck effect is addressed in the next 
section. 


A. Two-terminal setup 


Let us consider a two-terminal conductor, i.e., a system 
with two reservoirs [Fig. |2][a)]. This setting is relevant 
for spin-current rectification effect that emerges due to 
nonlinear spin transport in a conductor. Fig . Hljb) illus¬ 
trates the spin-current rectifier from Ref. j^. It includes 
a quantum wire with an asymmetric potential barrier in 
a uniform magnetic field. 

First of all, we notice that charge conservation implies 
AQi -I- AQ 2 = 0. This leads to /i -I- /2 = 0. Therefore, 
we have three independent currents, h, Isi, Is 2 - (Recall 
that spin may not be conserved in the conductor. Isi 
and Is 2 are defined at the interfaces between the con¬ 
ductor and terminals.) Besides, gauge invariance tells 
us that transport quantities depend only on the volt¬ 
age difference Vi — V 2 between the reservoirs and not on 
the absolute values of their chemical potentials. Thus, 
without loss of generality, we are allowed to set V 2 = 0. 
Therefore, we have three independent biases Fi, Fsi, Vs 2 
to drive the currents. 

Consider linear transport for a warm-up exercise. 
From the general fluctuation relation (l29ll . we derive in 
Appendix ID II the following relations for the noises 


5n(B) = 4T 
5'fi(B) = 4T 


5/1 (B) 
dVi ’ 
dlsijB) 
dVsi 


^“(B) = 2T 


dhjB) 

dVsi 


dlsijB) 

dVi 


(30) 

(31) 

(32) 




quantum wire 




H 



FIG. 2: (a) Two-terminal setup for charge and spin trans¬ 
port. The arrows represent our convention about the positive 
directions of the currents, (b) A spin-current rectifier [from 

Ref. [ 13 . 


and the linear conductances 


5/1 (B) 
dVi 

dh{-B) 
dVi ’ 

(33) 

dIsi{B) 

dVsi 

5/si(-B) 
dVsi ’ 

(34) 

dh{B) 

a/si(-B) 

(35) 

dVsi 

dVi ’ 


where all quantities are evaluated at equilibrium V = 
V 5 = 0. The relation (1501) is the famous Nyquist for¬ 
mula, and (I31II32I) are its generalizations for spin currents. 
The relation (l33l) is the Onsager reciprocity relation, and 
(I34II35I) are its generalizations for spin currents. Of par¬ 
ticular interest is the relation (I35|) . It has an additional 
minus sign compared to (1551) and (1551) . This additional 
sign originates from the fact that spin is odd under time 
reversal. 

Let us now turn to the nonlinear transport. We show 
in Appendix ID II that 


7lll(B) : 

"ni(B) : 

Aii(B) : 
"fii(B) : 


a^iAB) 2 92/1 (B) 

o^ g^fl(B) _ ,rr.2 d^Isi{B) 
dVsi ’ 

dVi ’ 

QC-S,-|- 

- 

dVsi 


and 


dS+ _ 

dVi 


ac"S'.+ 

dVsi 


dVl ’ 

5F|i ^ 
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dVi 


= 6T 


, 92 /- 


dVsi 


dVl ’ 


(36) 

(37) 

(38) 

(39) 

(40) 

(41) 


where we use the notation = 0(B) ± 0(— B) for 
convenience. Furthermore, we show in Appendix lD ll that 


_a2/i(B) dS,^{-B) ,^d^h{-B) 

dv^ dVi dv^ ’ 

d^IsijB) ^ dSf,{-B) _ dHsi{-B) 

5F|i dVsi dV§-^ 


_ = ^>gn(-B) _ d^Isi{-B) 

dVgi dVi dVidVsi 

d^IsijB) ^ dSn{-B) _ dH,{-B) 

dVi dVsi dVsidVi 


(42) 

(43) 


(44) 

(45) 


Again, all quantities are evaluated at V = Vs = 0. One 
can derive an infinite number of additional fluctuation 
relations. We focus on the above results because of their 
simplicity. In particular, only currents and second noise 
cumulants enter Eqs. (I40II45I) 




















































The relations (l36ll . (l38l) and (1401) for the third cumulant 
and nonlinear transport coefficients have been derived 
beforeii, and the relations §3), (IMl) and (mu are their 
generalizations for spin currents. Note that according to 
(l38l) . (7111(6) is odd under B —^ —B, while (7fii(B) is 
even. We notice that the right hand sides of Eqs. (I42II45I) 
can be obtained by differentiating the differences of the 
left and right hand sides of the Nyquist formulas (I30I31D . 
Therefore, Eqs. (I42M5I) can be understood as relations 
between the nonlinear conductances and the degree of 
the violation of the Nyquist formulas away from thermal 
equilibrium. Finally, it is worth to mention that the odd 
parity of spin under time reversal leads to various differ¬ 
ences (not just different signs on the left!) between (|38l) 
and (l3^ . (l40ll and dS]), (l42|l and (|4^ . and (l4^ and (|45|l . 


B. Multi-terminal setup 


The fluctuation relations for two-terminal setups pre¬ 
sented in the previous subsection can be generalized to 
multi-terminal setups. Below, we list the relations anal¬ 
ogous to the above two-terminal relations (see Appendix 
ID21 for derivations). The three-terminal geometry, il¬ 
lustrated in Fig. [21 is particularly interesting due to its 
relevance to transistors^2i2S. 

First of all, the generalizations of the Nyquist formula 
read 


5y(B) 

^y(B) 

5“(B) 


2T 

2T 


. dV, 
r5/s.(B) 
. dV, 


2T 


r5A(B) 

. dVsj 


diA^y 

dV, \ ’ 

dV, 

disA^y 

dV, 


(46) 

(47) 

(48) 



FIG. 3: Sketch of a three-terminal setup. 


The Onsager reciprocity relations read 


54(B) 

54 

54 (-B) 

54 ’ 

(49) 

54* (B) 
dVs, 

54 , (-B) 

5Vs* ’ 

(50) 

54(B) 

5Vs, 

54 , (-B) 

54 ■ 

(51) 


We emphasize the minus sign on the right hand side of 
dSlD. Relations, analogous to (I36I37|1 . read 


a,fc(B) =T 


dVk 


„ ^ ^ dSyjB) 

C^fc(B) =T^^-2T^ 


dVjdVk 


4- c.p; (52) 

dVsk “ dVsjdVsk ^ 

where “c.p.” stands for the cyclic permutations of the 
indices Relations analogous to (1551) and (1551) read 


Cyfc(B) = r- 


dVk 


= T- 


dS. 


jk 


514 


= T- 


.dS, 


ki 


dV, 


dS, 


s,+ 




= T- 


ds: 


s,-i- 


jk 


f)QS,+ 
_ rpOOki 


dVsr dVsj 

and relations analogous to (l40|) and (l4T1) read 

52 /+ 


dsy 

TnT + = 2^ 

aVk 

dST, 


dVk 


+ c.p. = 6T 


dS: 


S.+ 


dVsk 


dS: 


■s,- 


c.p. = 6T 


c.p. = 2T 


dVjdVk 

52 /- 

dVjdVk 

dVsjdVsk 


C.p. 


c.p. 


dVsk ' ' \dVsjdVsk 
The relations analogous to 
524(B) _ 1 55 y(-B) 


c.p. 


c.p. 


dV,dVj 


2T 514 


524 ,( 6 ) 1 55g-(-B) 

^Vs^^VsJ 2T dVsk 


read 
, 524(-B) 

51454 

524 . (-B) 
dVskdVsj ' 


(54) 

(55) 

(56) 

(57) 

(58) 

(59) 


5%(-B) 
514514 ’ 
(60) 

524 , (-B) 
dVskdVsi 
(61) 


and the relations analogous to (15511551) read 
524 (B) _ 1 54^-(-B) 


2T 514 


^Vs^^VsJ 

d^IskjB) _ _ 

514514 2T dVsk 


524 .(-B) ^ 5244 -B) 


dVkdVsi 

(62) 


1 54 , (-B) 


dVkdVsj 

524 (-B) 524 (-B) 

514,54 514,514 


(63) 

All the quantities above are evaluated at V = Vs = 0. 
The fluctuation relations hold for any indices = 

1,..., V, where N is the number of the terminals {N = 3 
in Fig. m. 
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C. Currents of spin-up and -down electrons 


We have focused above on the spin and charge currents 
and noises. It is also of interest to consider the currents 
and noises of spin-up and -down electrons. The latter 
currents are just linear combinations of electric and spin 
currents. We derive some fluctuation relations in this 
language which facilitates the comparison of our results 
to those of Ref. [13 in the next subsection. 

We will denote the currents of spin-up and -down elec¬ 
trons as lien where a = -I- corresponds to the spin-up 
current and a = — corresponds to the spin-down cur¬ 
rent. The charges of the spin-up and -down electrons in 
reservoir i are Qia = Qij2 — aeSzi. The conjugate chem¬ 
ical potentials Via = Vi — ^Vsi- Under time-reversal, 
Qia ^ Qicf! UcK ^ Ua; where -j- — , — -p. The 

fluctuation theorem (fT^ becomes 


P(AQ,V,B) =exp 


V AQ' 

T 


P(-AQ,V,-B), 


(64) 

where Q and V stay for the vectors with the compo¬ 
nents Qia and Via, Q and V are the vectors with the 
components Qia and Via, and V • AQ = J^ia ^ia^Qia- 
Below we derive two relations that connect the nonlin¬ 
ear conductance with the noise and the third cumulant. 
We sum the right and left hand sides of Eq. (IMl) over all 
possible AQ. We get (1)v.b = 1 on the left and expand 
the sum on the right to the third order in V. We find 


Cia,jp,k'i (V = 0, B) = T 

2 y 2 


95,/3,fc^(V = 0,B) 
^V^a 

9^4^(V = 0,B) 

dViadVjp 


f c.p. 

c.p. 


(65) 


where c.p. stays for the cyclic permutations of the pairs 
of the indices ia, j/S and kj. 

We next use Eq. (IMll to compute (AQia)v,B — 
{AQia)^ _-Q. We then subtract the right hand side of 
the equation from its left hand side and expand to the 
second order in V. The result is 


= 0,B) - aa,jPM(y = 0, -B) = 

^r55,a,,;3-(V = 0,B) 54a,fe4V = 0,B) 


dVk^ 

dS,a,jp{V = 0,-B) 


dV^-p 
= 0,-B) 


dVi 


k~f 


dV, 


iP 


( 66 ) 


D. Comparison with existing results 

In this section, we compare our result with the existing 
work. Refs. Hi] and [13, on fluctuation relations for spin 
currents and other multi-component currents. 

Ref. [U investigates a multi-component current of time- 
reversal-invariant quantities. At the same time, we em¬ 
ploy the fact that spin is odd under time reversal. This 


fact manifests itself in the sign factor (—I)^i("i+^») in 
the general expression (l29)) and in the sign factors in 
the fluctuation relations (1551) . (1511) . etc. If spin were even 
the sign factor in (1551) would be (—Note 
that Ref. [13 derives a similar relation to (lD7l) but with 
no minus sign [the minus sign in (ID7I) is carried over to 
(1551) ]. This illustrates the difference of the fluctuation 
relations for spin currents from the fluctuation relations 
for multi-components currents of time-reversal-invariant 
quantities. 

In Ref. [13, the authors derived several fluctuation re¬ 
lations for spin current without using microreversibility. 
In contrast, our derivation depends heavily on microre¬ 
versibility. For a detailed discussion of microreversibility, 
see Appendix [A] 

Our relatio n (1651) is similar to the fluctuation relations 
from Refs. [13S 3 and could be obtained without the use 
of microreversibility. On the other hand, Eq. (1551) is 
new and could not be obtained with the methods of Refs. 
I 22 II 23 I . Indeed, summing Eq. (1551) over all possible choices 
of the spin indices and taking the limit V = V, one 
reproduces a fluctuation relation in Ref[I3 for electric 
currents. The latter relation cannot be obtained without 
microreversibility and the same is true for our result (1661) . 


V. THERMOSPIN TRANSPORT 


We have assumed identical temperatures in the reser¬ 
voirs in the above discussion. In this section, we consider 
the case of the reservoirs at different temperatures. The 
temperature gradient leads to the spin Seebeck effect^i^— . 
In the ordinary Seebeck effect, an electric current flows 
between two conductors at the same voltage but different 
temperatures. In the absence of the symmetry between 
the two projections of spin, the currents of spin-up and 
-down electrons in response to a thermal gradient are not 
the same. Thus, a spin current is generated. 

Let us start with extending the fluctuation theorem 
(IT 2 I) to the case of the reservoirs with different tempera¬ 
tures. Consider the same setup as in Sec. m a conductor 
attached to several large reservoirs. The reservoirs are 
maintained at different temperatures Ti, different electro¬ 
chemical potentials Vi, and different spin chemical poten¬ 
tials Vsi. We follow the same measurement protocol as 
before: we measure the total energy, total charge and 
the z-component of the total spin of each reservoir be¬ 
fore and after an evolution period r. Following a similar 
derivation to the one in Sec. m one can show that 


P(AE,AQ,AS„T,V,Vs,B) 
P(-AE, -AQ, AS,, T, V, -Yg, -B) 


= exp 


'y ] {PiAEi -\- ^iAQi + (^iASzi) , 


(67) 


where P is the probability distribution function for ob¬ 
serving the change AEi of the energy, the change AQi of 
the charge, and the change ASzi of the z-component of 























10 


the spin in reservoir i during the evolution period r, with 
given Ti,Vi,Vsi and the magnetic field B. The quanti¬ 
ties are defined as Pi = 1/Ti, ^i = —VijTi and 

Q = -Vsjn 

We are now in the position to derive fluctuation rela¬ 
tions for thermospin transport. One can derive as many 
relations as in Sec. m We will focus only on the rela¬ 
tions that crucially depend on the odd parity of spin. As 
a warming up exercise, we consider the Onsager relations 
in linear transport. A detailed discussion of the Onsager 
relations in the spin Seebeck effect can be found in Ref. 
0- The approach of Ref. M is very different from ours 
and builds on the Landauer-Buttiker formalism. 

The average heat current, injected into the conductor 
from reservoir i, is given by 

Ihi = - lim -{{AEi — ViAQi — VsiASzi))B,T,v.Vs■ 

T—^OC T 

( 68 ) 

This definition reflects the fact that only the excess en¬ 
ergy change AEi — ViAQi — VsiASzi is dissipated into 
heat; ViAQi + VsiASzi can be thought of as the potential 
energy. We also find it convenient to define the energy 
currents 

Isi = — lim — (A£’i)B,T,v.Vs- (69) 

r—^■oo T 


Using the definition of Ihi and the connection between 
Pi, (i and Ti, Vsi, we obtain the relation 

1 d^hkjB) ^ 1 a^g-(-B) 

T^Vs^^VsJ 2T dTk dTkdVsj ^Tk^Vs^ 

(73) 

where the derivatives are computed at Vi = Vsi = 0 
and Ti = T. During the derivation, several terms are 
cancelled due to the Onsager relation (1341) . By similar 
arguments, it follows from (1551) and (H51) that 

^92/sfe(B) , 9/sfc(B) 

dT,dT, + dT, 

1 9^,^,-(-B) d^IhP-^) d^hA-^) 

2T2 dVsk dVskdTj dVskdT ^ ’ 

where is the cross noise between the heat currents 
Ihi and Ihj- The right hand sides of (I73p and (|74l) differ 
by a minus sign, which again manifests the odd parity 
of spin under time reversal. We note that there is an 
additional term proportional to Sij on the left hand side 
of (1751) compared to dZSl). 

VI. SUMMARY 


Let us now consider Eq. dnzi) in the limit of zero voltage 
bias, ^i = 0. Eq. (l67l) has now the same structure as Eq. 
m with Pi in place of —Vi/T and 0 in place of —VsijT. 
It follows from the analogy with Eq. (l35|) that 

diEm _ dlsA-^) 

dQ dp, ^ ^ ’ 

where the derivatives are taken at equilibrium: Vsi = 0 
and Ti = T. Making use of the connection between Pi, Q 
and Ti,Vsi, the definition of Ihi and the fact that spin 
currents vanish in equilibrium, we obtain the Onsager 
relation 

dhm 9/s,-(-B) 

dVsj dT, ■ ^ ’ 

Note that the right hand side describes the linear spin 
Seebeck effect, while the left hand side describes the lin¬ 
ear inverse spin Seebeck effect. 

We now turn to the nonlinear Seebeck effect and thus 
go beyond Ref. M Below we establish relations analo¬ 
gous to dMl) and dMl). Again, we make use of the simi¬ 
larity between the fluctuation theorems (lT2l) and (l67)) in 
the case of = 0. Taking Eq. (I62p and making the sub¬ 
stitutions —Vi/T —>■ Pi, —Vsi/T (i and A —>■ lEi, we 
find that 

d^lEkiB) i55,§(-B) d^Is,i-B) d^Isji-B) 

dQdQ 2 dPk dPkdQ dPkdC, ' 

(72) 


In conclusion, we prove a general fluctuation theorem 
for spin currents. It imposes a great number of restric¬ 
tions on transport coefficients and fluctuations of the 
spin, charge and heat currents in linear and nonlinear 
transport. We have focused on the fluctuation relations 
for currents, noises and third cumulants, and their lower 
order derivatives. Several relations for the transport, 
driven by the gradients of chemical potentials, are derived 
in Sec. II V Al and Sec. IIVBI Sec. IVl addresses thermospin 
transport. 

Our derivation relies on the combination of microre¬ 
versibility with the assumption that the ^-component of 
spin conserves approximately in large parts of the reser¬ 
voirs near the conductor. If the latter assumption does 
not hold then it may not be even meaningful to speak 
about the spin currents, injected from the reservoirs into 
the conductor. Note that very long spin-coherence times 
have been reported for donor spins in silicon^. 

Our results do not depend on a particular model and 
thus apply to many systems. One example is the spin 
current rectifier— Since our fluctuation relations are 
exact, they can be employed for testing approximations 
in theoretical calculations. 

The discussion in the preceding sections is most di¬ 
rectly connected with mesoscopic conductors that can 
carry both spin and electric currents. At the same time, 
the presence of an electric current is not essential for 
the validity of our results. They also apply to pure spin 
currents. Eor example, our fluctuation relations hold for 
spin diffusion in insulators^. The above calculations as¬ 
sume that the spin current is carried by particles with 
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spin 1/2. It is straightforward to extend our results to 
higher spins. Another interesting setting is topological 
matter with chiral transport^. We expect our fluctu¬ 
ation relations to simplify considerably in such systems 
since chirality implies the vanishing of many transport 
coefficients^^. A simple example is discussed in Appendix 
11 
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Appendix A: Microreversibility 

The principle of microreversibility is expressed by Eq. 
(El)- Its physical meaning is simple: after reversing the 
directions of the magnetic field and the velocities and 
spins of all particles, the system traces its evolution back¬ 
wards in time. This is a straightforward consequence of 
the laws of quantum mechanics. At the same time, the 
calculations2^i2i, based on the Landauer-Biittiker formal¬ 
ism, seem to conflict with the principle. In this section 
we clarify the origin of the conflict. We will see that 
while the Landauer-Biittiker formalism works well in the 
linear response regime, its application to fluctuations in 
nonlinear transport faces challenges. 

The Landauer-Biittiker formalism^ uses a single¬ 
particle language. Electrons are treated as non¬ 
interacting particles in a self-consistent electrostatic po¬ 
tential. The self-consistent potential is different for the 
opposite directions of the magnetic field^i^. The ori¬ 
gin of that asymmetry is illustrated in Fig. [4] In that 
example, charged particles enter the conductor through 
terminal 1. Their trajectories depend on the magnetic 
field. Depending on its sign, the particles exit through 
terminal 2 or 3. The average charge density p(r,B) 
is nonzero only along the particle trajectories. Hence, 
p(r,B) ^ p(r,-B). Since the self-consistent electrostatic 
potential (j) is determined by the charge density, 

<^(r,B)^</(r,-B). (Al) 




FIG. 4: Trajectories of particles in opposite magnetic fields. 


The central quantity in the Landauer-Biittiker formal¬ 
ism is the single-particle scattering matrijs^. Refs. MM 
compare the scattering amplitude from the initial state 
with the momentum k,; to the final state with the mo¬ 
mentum k f for an electron in the magnetic field B in the 
self-consistent potential (/(r, B) with the “time-reversed” 
scattering amplitude from the state with the momen¬ 
tum —k/ to the momentum —in the field —B in the 
potential (/(r, —B). If the two self-consistent potentials 
B) and (/(r, —B) equaled the two scattering ampli¬ 
tudes would be complex conjugate. Eq. (EH) implies 
that they are not in fact complex conjugate. This was 
interpreted as microreversibility breaking^^. Indeed, the 
calculations^i^i with such scattering amplitudes contra¬ 
dict the results from microreversibilityi^. 

The origin of this apparent breakdown can be best un¬ 
derstood if we forget about quantum mechanics and con¬ 
sider a classical many-body system. Let us first use a self- 
consistent approximation. We find a self-consistent elec¬ 
trostatic potential (j){r, B) and consider separately the 
motion of each electron in this potential. Let us now 
reverse the direction of the final velocity of the elec¬ 
tron and the direction of the magnetic field B —>■ —B. 
If the self-consistent potential remained the same then 
the electron would trace the same trajectory backwards. 
A different self-consistent potential ^(r, —B) implies a 
different trajectory, i.e., microreversibility breaks down. 
How does this compare with the exact solution of the 
equations of motion of a many body system? The po¬ 
tential energy of a system of electrons depends only on 
their positions and does not depend on the magnetic 
field. It follows then easily that after one reverses the 
magnetic field and the velocities of all electrons then 
each electron traces its trajectory backwards, exactly 
as the microreversibility principle demands. In other 
words, microreversibility breaking is an artifact of the 
single-particle self-consistent approximation inherent to 
the Landauer-Biittiker formalism. Microscopic calcula¬ 
tions for many-body quantum-mechanical models beyond 
the Landauer-Biittiker approximation^Ii^^ do agree with 
the predictions^^ from microreversibility for charge trans¬ 
port. The goal of this work is to investigate the conse¬ 
quences of microreversibility for nonlinear spin transport. 

The above discussion shows that a magentic field does 
not destroy microrevesibilty if it is present without a 
magnetic field. Microrevesibility without a magnetic field 
has been addressed by many author s The results, 
based on microreversibility, can be expected to apply in 
a Hamiltonian system, where the conductor can only ex¬ 
change energy with the reservoirs. This condition cannot 
be relaxed if the reservoirs are maintained at different 
temperatures. When all temperatures are the same such 
energy conservation is no longer necessary. Indeed, the 
insulating enviroment can be included into one of the 
reservoirs in all derivations. Another condition remains 
crucial for the validity of our results even when all tem¬ 
peratures equal as discussed in Refs. US the charge 
transfer must be small in the detector that probes the 










12 


(a) (b) 



FIG. 5: Detecting spin chemical potentials, (a) A qnantnm 
wire proximate to the reservoir as a probe, (b) Schematics 
of the energy snbbands of the wire and the relative values of 
different energies. 

system. Otherwise the fluctuation relation is modifed^. 


Appendix B: Detection of spin currents Is and 
chemical potentials Vs. 

In order to connect our results with spintronic exper¬ 
iments, we briefly review the ways to measure the spin 
currents Isi and spin chemical potentials Vsi in this ap¬ 
pendix. 

Detecting spin currents is of crucial importance for 
spintronics and numerous methods to do so where im¬ 
plemented and/or proposed. Examples include the 
use of the spin-current induced Hall effect^, magnetic 
resonance^ and optical techniques^. The measurement 
of the chemical potentials Vs reduces to the measurement 
of electric and spin currents. Below we discuss one way 
to extract Vs from currents. 

The setup is illustrated in Fig. [5l A quantum wire is 
brought close to reservoir i so that electrons can tunnel 
between the wire and reservoir. As discussed below, the 
tunneling current contains information about the chemi¬ 
cal potentials = —eV-f and /i^ = —eVj, of the spin-up 
and -down electrons in the reservoir. Vs can be computed 
from Vs = - in,. 

Since we are interested in a situation with an external 
magnetic field, the lowest subband in the wire is split 
into the spin-up and -down subbands [Fig. [Sj^b)]. Let us 
assume that the band bottom Ff for spin-up electrons is 
lower than the band bottom for spin-down electrons. 
We will assume that the electro-chemical potential p, of 
the wire is the same for the spin-up and -down electrons 
and stays between E-f and E^. We will also assume that 
the temperature T <C Ei — Ef. Then only the lower 
spin-up subband is populated. Our main focus is on the 
limit of small Vs. By changing the electric potential of 
the wire, the chemical potential /i can be made close to 
/i|. In view of the smallness of Vs this implies that the 
chemical potential ^ 4 , of the spin-down electrons in the 
reservoir is also close to ii. In such situation we expect 
no tunneling of spin-down electrons between the wire and 
the reservoir. On the other hand, the tunneling of spin- 
up electrons disappears only at /r = /i|. Thus, a zero 


total tunneling current indicates /i = /Xf. This gives a 
way to measure the chemical potential /i| solely from a 
measurement of an electric current. 

In order to find Vs one also needs to measure ^ 4 ,. This 
can be accomplished in one of two ways. One can use 
another quantum wire with an opposite ( 7 -factor. In such 
wire the order of the spin-split subbands reverses and 
thus an electrical current measurement yields Alter¬ 
natively, one can work with the same quantum wire as in 
the hrst experiment. By changing the electric potential 
and the charge density of the wire, one drives the chem¬ 
ical potential 7 / so that ii lies above both band bottoms 
E^ and E^. In such situation and focusing on the regime 
fi fa Ki there is a tunneling current of spin-down 
electrons unless ii = In that point, the tunneling cur¬ 
rent is fully spin-polarized. The polarization of the carri¬ 
ers can be determined by a simultaneous measurement of 
the spin and electric currents between the reservoir and 
wire. This can be used to deduce 7 / 4 ,. 

Appendix C: Higher-order cumulants and 
generating function 

In this appendix, we briefly describe higher-order cu¬ 
mulants of the heat, electric and spin currents. We con¬ 
sider a general case that the reservoirs have different tem¬ 
peratures, as well as different electric and spin chemical 
potentials. In general, cumulants can be defined through 
the generating function, given by 

E{x,y,z,T,V,Vs,B) = 

lim - In I V g-(E> 

r^oo T 1 ^ 

L AE,AQ,AS2 

xP(AE,AQ,AS„T,V,Vs,B)|. (Cl) 

where the bold symbols x,y, z,T,V,'Vs are vectors, 
e.g., X = {xi,X 2 , ■ ■ ■ ,xn), and N is the total number 
of the reservoirs. 

The nth order cumulants of heat currents are given by 

Ct.k = id,. - VA, - VsA.m, - V,dy, - Vs,a.,) 

.. .{d,, -Vkdy, -Vskd,,)E(x,y,z,T,V,Vs,B) 

(C2) 

where x,y,z are eventually set to 0 after the derivatives 
are taken. The number of the indices i, j,..., fe is n, and 
the indices may or may not be the same. The nth order 
cumulants of electric currents are given by 

= dy,dy, . . . dy,T{X, ?/, 2, T, V, Vg, B) (C3) 

and the nth order cumulants of spin currents are given 
by 

y, z, T, V, Vs, B) (C4) 
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where again x,y,z are eventually set to 0 after the 
derivatives are taken. In general, one can define various 
cross cumulants of heat, electric and spin currents by cer¬ 
tain combinations of the derivatives {dxi—Vi dy^ — Vsi ), 
dyj and 9^^,. Note that the noises and third cumulants 
in (fTTra differ from the above definitions by a conven¬ 
tional factor of 2. 

General fluctuation relations can be obtained for 
higher order cumulants, though they are less relevant to 
experiments. The fluctuation theorem (1671) can be trans¬ 
lated into a symmetry of the generating function, given 
by 


where /3 is the vector with the entry [3i = \/Ti^ ^ is the 
vector with the entry = —Vi/Ti^ and C is the vector 
with the entry Q = —YsilTi- To obtain the fluctuation 
relations, one can perform the Taylor expansion of the 
generation function T on both sides of Eq. (IC5D around 
X = y = z = Y = V5 = 0 and Ti = T. Comparing 
the Taylor coefficients on the two sides, in principle, one 
can obtain an infinite number of fluctuation relations. 
In the main text, we instead use Eq. (1291) to derive the 
fluctuation relations for the lower order cumulants. 


T{x,y,z,T,V,Vs,B) = 

T{l3-x,^-y,z-C,T,V,-Vs,-B), (C5) 


I 

Appendix D: Derivation of fluctuation relations from Section IV 


In this appendix, we provide derivations for some of the relations, listed in Sec. CYl from the general expression 

(EH). 


1. Derivation of two-terminal relations 


To prove the two-terminal relations (13111) - (l45|l . we start with a simple helpful identity (AQ°AS°) = 1. As a 
consequence. 


= Su.oSfj^o- (Dl) 

As mentioned in Sec. IIV A[ in the two-terminal setup, charge conservation makes AQ 2 dependent on AQi and the 
gauge invariance principle allows us to set V2 = 0. So, we can omit the indices n2 and 1 ^ 2 , and use a shorter notation 
instead of the full notation L"’™. 

Let us derive (l30l) and (|^ . Taking the left hand side of (l29)l to be Lqqq(B), L}Qg(B), L2oo(^) respectively, we have 

Lg™(B)=L2°°(-B), (D2) 

L}™(B) = -Ll™(-B)-Lg™(-B), 

^2oo(B) = L^2ooi-^) + i}oo(-B) + IlI°°{-B). 


Making use of dnn), we obtain 


l200(B) = -2L}°0(B), (D3) 

L™(B) = Li™(-B). (D4) 


With the definition (l28l) of L and the definitions given in Sec. IIII A1 we can identify Lqqo = <S'ii/2, Lj™ = —Tdh/dVi. 
One immediately obtains (l30l) . (|33p from (|D3p . (jD4l) . Similarly, we obtain (iMl) . (IM)l by taking the left hand side of 
(EH) to be L020 (b),loio(B),LO™(B). 

To derive (|^ and (IMI) . we take the left hand side of (l29)l to be LoJq(B), ^^^(B), L?og(B) and L5 io(B) respectively. 
We then have 


Lll°{B) = -Lll°i-B), 

LJ?°(B)=LJ?°(-B) + L“0(-B), 

l?J°(b) = l;j°(-b) + l“°(-b), 

0 = L??°(B) = -T??°(-B) - Li™(-B) - L?J°(-B) - L“°(-B), 


(D5) 








and hence 
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L“°(B) = -LJ?°(B)-L?J°(B). (D6) 

L™(B) = -L?J°(-B). (D7) 

With the identifications LJqq = Sfi/2, = —Tdli/dVsi, and = —Tdlsi/dVi, we immediately obtain (15^ . 
from dnSl), (iDTll . 

As an aside, we find that S'ii(B) = 5'ii(—B) and 5'fJ(B) = —B) from the equations (ID2II and (ID5D respec¬ 
tively. A similar relation 5f2(B) = S'f2(—B) holds. 

Next, we derive relations for nonlinear coefficients and the third cumulants. We first prove (l36l) . (l38l) . (ITO) and 
(1^ . Taking the left hand side of (1^ to be Looo(B)! ^ioo(B), A20o(^)i-^30o(^) respectively, we have 

rSOO/'-D'v _ 7-300/ TjN 

-^oool-’^'j ~ -^oool 

L?™(B) = L?™(-B) + l300(_b), (D8 ) 

Li™(B) = -LZ(-B) - L?°°(-B) - 
0 = L°™(B) = L°™(-B) + L^™(-B) + iL?™(-B) + iAg™(-B), 


With the equation m and with some computation, we arrive at 

l300(B) = -3A?™(B)-6L™(B), 
L?°°(B) + A?™(-B) = -2 [Li°°(B) + ^^(-B)] , 
L?o°(B) - L?™(-B) = -6 [Lio°(B) - l1™(-B)] , 
2Li™(B)=L?°°(-B)+4Li™(-B). 


Recognizing that 2AQQg = —Cm, 2L^gQ = TdSn/dVi, and 2L200 = —T'^d'^Ii/dVi, we immediately obtain (IM)) . ([38]), 
(HUl) and from (IDSp and the above four equations. 

To derive (lT7|l . (l39l) . (ITO and (l43|) . we consider the following consequences of (1^ 

7-030/t:>\ _ 7-030/ -ON 

-^0001-*=*^ ~ -^oool 

7-020/t:>\ _ 7-020/ idn 7-030/ td^ 

-^oiol-*^J ~ -^oiol -^oool 

io2o(B) = + ^i^ooo(-B), 

7-000/t:>\ _ 7-000/ xy\ rOlO/ r:>\ ^ 7-020/ r:i\ ^ r030/ t::>\ 

-^030V-*=»J — -^030V ^ 020 \ 2 ^ 010 \ q ^ 000 \ 


Similarly, we are able to obtain 


L°20(B) + L°?°(-B) = -6 [L°i°(B) + L°i°(-B)] , 
L020(B) - A°?°(-B) = -2 [l 010(B) - T°i°(-B)] , 
2L°i°(B) = -i°™(-B)-4L°i°(-B). 


Then, the relations dSTj), (|39l) . (ITO and (l4^ immediately follow. 

Finally, we prove (11^ and (H51) . Consider the consequences from (1^ 

4™(B) = -Aj™(-B) - L“0(-B) - 

i°°°(B) = L™°(-B) + l;]°(-B) + iL?g°(-B) + Li°°(-B) + L“0(-B) + iLj2°(-B), 
L°J0(B) = L«J°(-B) + L“°(-B) + 

i0™(B) = -L°;°(-B) - Ll%{-B) - hl%{-B) - L°J°(-B) - L“°(-B) - iLgJO(-B). 
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It follows that 


(D9) 

LOio(B) = -Lioo(-B) - iL200(-B). (DIO) 


Recalling that 2Ll°° = L^lo = -T^d^Isi/dVidVsi, 2 L ?20 = TdSf,/dVi, 2L%°o = -T^d^IsildV,^ 

L)°q = —T'^dli/dVidVsi and 2Lg°g = TdSn/dVsi, we immediately obtain (11^ and (H51) . 


2. Derivation of multi-terminal relations 


The derivations of the multi-terminal relations (I46l) - (l63p are very similar to the derivations in the two-terminal case. 
Below we show derivations for some representative relations. The rest can be derived in essentially the same way. 

Let us derive the relations ((4^ and dSU). Let i be an integer vector with the z-th entry being 1 and the other entries 
being 0. The notations j, fc,... are similarly defined. According to the general expression ([29|), we have 

= i-i(-B) + (Dll) 

4° (B) = 4(°(-B) + 4’;i,(-B), (D12) 

L°;(B) = -L°;°(-B) - 4’°(-B) - L%(-B) - 4(i(-B). (D13) 

Using the identity m, we find 

4(^(B) = -4(°.(B)-L°^(B), 

4“(B) = -L°;^'(-B). (D14) 

Recall that Lg-j, = S^/2, Lq°. = —Tdh/dVsj and L°q = —TdIsj/dVi. We immediately obtain (l48)) and (IbTl) . 

Let us also derive dSSl), dSll), (Ull), (IMl) and m- For simplicity, we assume i ^ j ^ k. However, the final result 
does not depend on this assumption. According to the general expression (IMl) . we have 


l: 

L] 


0,i+j 

O.k 

0,i 

OJ+k 


tO,0 

^0,i+j+k 


(B) = 

(B) = L°:;+fc(-B) + L°;t+^(-B) + Llf^i-B) + L°f^+^{-B), 
(B) = .^,(-B) - L°;}+,(-B) - L°:i+,(-B) - L°;t,,(-B) 


T 0.2+J 

^O.fc 


(-B)-L 


0,j+fe 

0 ,i 


(-B)-L 


0.2+fc 

0,J 


(-B)-L 


0,0 


(-B). 


(D15) 


Using the relations, obtained from the above by cyclic permutations of the indices i,j, k and the identity (IDID . it is 
not hard to show that 


rO.i+j+k/ 


.0,0 =(B) = -L°;;.^,(B)-L";r^( 




B) 


c.p., 


where c.p. stands for cyclic permutations of the indices i, j, fc, and 


■^(B)+c.p. 

L°;t+^(B)+c.p. 


L 


o:i""^'(-B) 


L 


0,2+J 

O.fc 


(-B) 


C.p. 

c.p. 


3 { [^o;^-+JB) + C.p] + J-B) + C.p] } 


{ 


L 


OA 


L 


O.i 


0,j+k 

_ rOJ+k 


(B) 


c.p. 

rO,j 


T 0 ,i 
^0,j^k 


(-B) -(- c.p. I 


oy+JB) = - L^;r"(-B) - 


Recall that 2Loo^^'^^ = 2Lo'^’^^ = TdSf^/dVsk and Lq = —T^d'^Isi/dVsidVsj. We immediately obtain 

(l53ll . (l58l) . (l5^ and (IHTI) . The relation (1551) follows from (ID15I) . 

Finally, we derive the relation (IMl) . According to the general expression (IMl) . we have 


k,0 

^O.i+j 

0.0 


(B) = ,(-B) - L^'](-B) - LSj(-B) - LS.r^(-B), 


- kA 


k;j ( 


. k.i+j , 


^w(B) = L^“+,(-B) + Afc;}(-B) + L"’j(-B) + L"’;o^^(-B) 


0,i 


0,3, 


-OA+j, 


+ 4f+i(-B) + AS;}(-B) + ASj(-B) + LS,r^(-B). 


- kA 


k,j , 


- kA+j ( 


(D16) 

(D17) 
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where we have assumed i ^ j for simplicity, but the final result holds even ii i = j. Using the identity m and 
adding the above two equations, we obtain 

.(B) = + L“;;.(-B) + L°’j(-B). (D18) 

Recalling that = TdSfj/dVk, Lo'°+j = -T’^d’^h/dVsidVsj and = -T'^d^Isj/dVkdVsi, we immediately 

obtain (l62l) . 


Appendix E: Spin transport in a quantum Hall 
system. 

In this appendix we briefly consider an example of a 
chiral system. In chiral systems there is transport in 
one direction only, e.g., clockwise. Chiral transport re¬ 
sults in a ‘stronger’ casualty then usual: Not only past 
does not depend on the future but also what happens up¬ 
stream does not depend on what happens downstream. 
Such stronger causality principle greatly simplifies the 
fluctuation relations and leads to stronger restrictions on 
transport coefficients than in non-chiral systemsi^Ti^S. 

Our example is illustrated in Fig. [6l We consider a 
two-dimensional electron gas in a strong magnetic field 
in the integer quantum Hall regime with the filling factor 
1 / = 1. The bulk of the sample is gapped and all transport 
occurs on the edges. There is also tunneling into conduc¬ 
tor C. The source S and the drain D are maintained at 
the same temperature T and voltage V. There is a volt¬ 
age and temperature gradient between S and C. We are 
interested in the spin current Ic, flowing into C, and the 
spin current noises Sd and Sc, detected in D and C. The 
current is fully polarized at zz = 1. Thus, the spin cur¬ 


rent Is = —Il2e. Hence, a non-equilibrium fluctuation- 
dissipation theorem for spin currents and noises can be 
deduced from the results for the noises of the charge 
currentsi^^— : 


Sd 


Sc + 


2T die 
e dV 


T 

2tt 


(El) 


Note that Eq. (EH) relates quantities defined at the same 
magnetic field. Note also that Eq. (EH holds in the 
nonlinear transport regime, i.e., at a finite voltage drop 

F-Uc- 


V 



FIG. 6: Three-terminal quantum Hall setup. 
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